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ABSTRACT 

This  is  an  analysis  of  three  dimensional  surface  waves  which  occur  in  a  vessel  of  rectangular 
cross  section  and  finite  depth.  The  linearized  problem  has  double  eigenvalues  for  particular 
combinations  of  the  parameters.  It  is  shown  that  eight  solution  branches  of  finite  amplitude  are 
emitted  by  the  double  eigenvalues.  Splitting  of  the  double  eigenvalues  resul  s  in  a  secondary 
bifurcation.  The  direction  of  the  emitted  branches  for  the  multiple  and  secondary  bifurcation 
changes  with  the  depth  of  the  fluid.  Finally  it  is  shown  that  the  formal  solutions  obtained  are 
not  uniformly  valid  and  an  additional  expansion  in  the  region  of  non-uniformity  shows  that  the 
wave  field  changes  type.  One  possibility  in  this  region  is  a  field  of  three-dimensional  cnoidal 
standing  waves. 
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SECONDARY  BIFURCATION  AND  CHANGE  OF  TYPE  FOR 
THREE  DIMENSIONAL  STANDING  WAVES  IN  SHALLOW  WATER 

Tlionihs  .1  Bridges 

1.  Introduction 

Three  dimensional  irrotational  standing  waves  of  finite  amplitude  in  a  basin  of  rectangu¬ 
lar  cross  section  and  finite  depth  are  considered  in  this  paper.  The  governing  equation  and 
boundary  conditions  in  dimensionless  form  for  the  potential  and  the  wave  height  are 

+  +  =  o  (i.i) 

di*  dy-^  ^  dz^  ^  ^ 


and  on  y  =  z,  t) 


=  0  on  solid  boundary 


drj  dt}  dd>  jdridd^  d4> 


(1.1) 

(1.2) 

=  0 

(1.3) 

-t- 1>  =  0 

(1.4) 

dt  dx  dx  dzdz  dy 

where  ^  =  ^,2a  is  the  vessel  length  in  the  x  direction,  2b  is  the  vessel  length  in  the  z  direction, 
^  ~  ia>  ^  i®  th®  vessel  still  water  depth,  f  =  x  'vhere  H  is  a  measure  of  the  wave  height. 

The  linearized  problem  for  (1.1)-(1.4)  has  eigenvalues  (bifurcation  points,  linear  natural 
frequencies) 

^0  ~  \/^^mntanhACTirt^  (L5) 


where 


Amn  =  JTv/m*  + 


The  solutions  for  finite  but  small  amplitude  emitted  by  the  simple  linear  eigenvalues  have 
been  found  by  Verma  iz  Keller  (1962)  using  a  formal  perturbation  expansion.  However  at 
particular  combinations  of  (r7i,n)  and  (  there  are  double  eigenvalues.  It  is  shown  in  section 
2  that  eight  solution  branches  are  emitted  by  these  double  eigenvalues.  This  is  done  using  a 
formal  perturbation  expansion  similar  to  Verma  k  Keller 
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In  an  interesting  discovery  Bauer,  Keller,  ir  Reiss  (1975)  observed  that  2is  a  multiple 
bifurcation  point  is  “split”  (by  varying  an  auxiliary  parameter  (in  this  case  0)  into  primary 
bifurcation  points  a  socofidary  bifurcation  may  txcur 

111  the  standing  wave  problem  there  are  four  <  branches  emitted  by  the  double  eigenvalues. 
Splitting  a  double  eigenvalue  results  in  two  simple  bifurcation  points  with  a  single  +£  branch 
emitted  from  each.  Rather  than  vanish,  the  other  two  branches  from  the  double  eigenvalue 
slowly  depart  by  creeping  up  a  primary  branch.  Bauer,  Keller,  ii  Reiss  developed  a  perturbation 
method  to  analyze  this  phenomena.  Subsequently  this  theory  has  been  successfully  applied 
by  Matkowsky,  Putnick  ii  Reiss  (1978)  to  the  buckling  of  rectangular  plates,  Kriegsmann  & 
Reiss  (1978)  to  the  theory  of  magnetohydrodynamic  equilibria,  and  has  been  extended  to  the 
bifurcation  from  triple  eigenvalues,  which  results  in  secondary  and  tertiary  bifurcation  by  Reiss 
(1983). 

This  theory  is  applied  in  Section  3  to  show  that  secondary  bifurcation  occurs  in  the  neigh¬ 
borhood  of  the  double  eigenvalues.  The  secondary  bifurcation  points  for  a  perturbed  square 
cross  section  are  found  as  a  function  of  6  and  the  mode  numbers  and  solutions  along  the 
secondary  branches  are  derived.  It  is  found  that  the  jump  to  a  secondary  branch  produces 
interesting  irregular  wave  forms. 

The  solution  of  Verma  ic  Keller,  the  finite  amplitude  solutions  found  in  section  2  emitted 
by  the  double  eigenvalues,  and  the  secondary  bifurcation  phenomena  elucidated  in  section  3 
are  valid  for  a  small  range  of  amplitude  only.  It  is  shown  that  when  <  =  0(8^)  the  higher  order 
terms  are  no  longer  of  higher  order  and  the  expansion  breaks  down. 

In  section  4  a  separate  analysis  is  performed  with  the  assumption  that  c  =  0(^^)-  As  a 
first  approximation  weakly  three  dimensional  waves  are  considered.  This  is  done  by  looking  in 
the  region  of  the  ^  -  6  -  (  parameter  space  where  the  triple  balance  holds 

-  0(e) 

^^^0(.) 

This  analysis  results  in  a  field  of  standing  K-P  waves.  A  set  of  two  non-interacting  (to  first 
order)  solutions  of  the  K-P  equation  (Kadomtsev-Pet  viashvili  1970).  The  K-P  equation,  which 
is  rich  in  solutions,  has  been  studied  in  some  detail  by  Dubrovin  (1981)  and  Segur  A  Finkel 


(1984).  Therefore  as  the  amplitude  increases  beyond  <  =  0(6*)  the  wave  field  changes  type. 
The  fact  that  there  are  multiple  solutions  when  r  <♦'  6*  and  the  fact  that  the  K-P  equation  is 
rich  in  solution  possibilities  suggests  that  tlie  region  for  <  -  0(6*)  will  result  in  a  “rats  nest 

of  branches  of  solutions 

The  results  with  the  K-P  equations  are  for  weakly  three  dimensional  waves.  A  further 
analysis  is  performed  in  section  4  with  this  tissumption  relaxed.  This  analysis  results  in  a  wave 
equation  to  leading  order.  The  solvability  condition  at  the  next  order  results  in  a  functional 
differential  equation  for  the  leading  order  term.  A  complete  solution  is  not  found  but  it  is 
shown  that  one  solution  is  a  set  of  four  oblique  travelling  cnoidal  waves  which  combine  to  form 
a  three  dimensional  standing  wave.  Distributions  of  the  wave  height  are  shown  for  these  waves. 
It  is  expected  that  this  equation  will  yield  other  interesting  possibilities. 

2.  Primary  Bifurcation  When  e  <  0(6*) 

Linearizing  the  set  of  equations  about  the  still  water  level  results  in  a  linear  problem  which 
has  eigenvalues  (linear  natural  frequencies) 


and  eigenfunctions 

Tji  =  cosQmicos/d„zsint  (2-2) 

,  _  1  cosh  A,nn(y  "^  ^)  —  n  —  ,  /o 

01  =  (To  ‘ - T-T - - —  cos  a„z  cos  cos  t  (2.3) 

cosh  A„„6 

where  A^n  =  v/oj,  +  f*/3*,  —  nn,  x  =  z+|,z  =  r-t-^,  and  m,n  are  the  mode 

numbers  in  the  x,z  directions. 

The  solutions  which  bifurcate  from  the  linear  eigenvalues  (2,1)  were  first  found  by  Verma 
Keller  (1962)  using  a  perturbation  expansion  in  the  amplitude.  They  found  that  the  natural 
frequency  for  the  three  dimensional  standing  wave  has  the  following  form  as  c  — *  0 


{3a*  +  Xl„  -  4al)^  {3a*o  +  . 

64<7^(a„,  tanh  2q„5  -  2ff^)  64a^{0„(  tanh  20„^6  -  2a^) 

Taking  q„,  -  -  I  and  ^  ~  i  this  agrees  with  the  expression  in  the  paper  by  Verma  ii. 

Keller.  The  bifurcation  for  the  frequency  is  sub-  or  supercritical  depending  on  the  value  of 
6.  As  6  ^  oo  (deep  water)  the  bifurcation  is  subcritical  and  as  6  — *  0  (shallow  water)  it 
is  supercritical.  Figure  2  of  the  paper  by  Verma  ii  Keller  shows,  for  the  first  mode,  that  as 
the  vessel  cross-section  departs  from  being  a  square  the  critical  depth  5’  (the  point  where  the 
bifurcation  changes  from  sub-  to  supercritical)  increases.  The  same  phenomena  occurs  for  the 
higher  modes  as  well  with  the  critical  depth  being  smaller  as  the  mode  number  increases  (for 
fixed  ^).  Consequently  when  a  lower  mode  is  supercritical  it  may  be  that  a  higher  mode,  with 
all  other  parameters  being  equal,  may  bifurcate  subcritically  suggesting  the  possibility  of  an 
intersection  of  the  branches  at  finite  amplitude. 

It  can  be  shown  however  that  the  asymptotic  solution  obtained  by  Verma  id  Keller  is  not 
uniformily  valid.  Taking  the  limit  as  6  —•  0  it  is  found  that 

02  k 

where  fc  is  a  constant  of  0(1).  Therefore  the  solution  is  valid  for  c  <  0(6*)  only.  In  section  4  the 
equations  will  be  reanalyzed  for  the  region  (  =  0(6*)  and  it  will  be  shown  that  the  solutions 
change  type  in  this  region. 

Inspection  of  (2.1)  also  shows  that  the  linearized  problem  has  double  eigenvalues  at  par¬ 
ticular  combinations  of  q^,  /?„,  and  For  example,  when  ^  =  1  every  pair  {am,0n)  such  that 
m  ^  n  is  a  double  eigenvalue  and,  for  example,  when  $  =  |  (m,n)  =  (1,4)  and  (2,2)  share  the 
same  eigenvalue.  In  fact  every  rational  $  will  have  an  infinite  set  of  double  eigenvalues.  It  will 
now  be  shown  that  these  double  eigenvalues  emit  multiple  branches  of  solutions.  For  brevity 
the  problem  of  a  square  vessel  ({  =  1)  will  be  considered  and  it  is  expected  that  the  analysis  at 
other  double  eigenvalues  will  result  in  a  similar  conslusion.  When  ^  =  1  the  bifurcation  points 
occur  at  (2  1)  with 

Amn  yam  (2  6) 

Therefore  any  pair  (m,n)  such  that  m  ^  n  will  result  in  a  double  eigenvalue.  The  analysis 
proceeds  by  formally  expanding  the  frequency,  potential,  and  wave  height  in  a  regular  pertur¬ 
bation  series.  The  leading  term  in  the  frequency  expansion  is  given  by  (2.1)  with  (2  6)  and  the 


leading  term  for  each  of  the  dependent  variables  is 


r)i  /4ii  jros5„;  -  ^i2COS,?„icosQm2isin< 


1  .  -  ^  ,  ,3-  .cosh  A„„(y  +  6) 

C)J  -  -  .■\iiCOSO„TCOS/<„C  .-ii-.ros/JnicosQ^z  - r— - - — cost 

ffii  coshA„„i5 

A  normalization  for  the  coefficients  is  chosen  such  that 


(2.7a) 


(2.76) 


A*i  +  A*2  =  1  (2.8) 

The  relative  magnitudes  of  Ah  and  A 12  are  determined  at  higher  order. 

Proceeding  in  the  usual  way  results  in  the  following  set  of  bifurcation  equations  after 
application  of  the  double  solvability  condition  at  the  third  order. 


[oiAh  +  a2AH  +  2a'2]Aii  —  0 
[a2AH  +  ajAj2  +  2ff2lAi2  =  0 


(2.9a) 


(2.96) 


which  along  with  (2.8)  form  a  set  of  three  equations  for  the  three  unknowns:  az,  An,  A12. 
The  coefficients  aj  and  aj  are  given  by 


gp  r  ,  „.(ai.  +  ^;;)  ,,2 


(3a^-3a»  +/9»)*  (3<73  -  3^4  + 


32ao[am  tanh  20^6  -  2<To]  Z2aQ\0n  tanh  2/9„6  -  2(7*] 


/?„)tanh(\/2(am  +  0n)6]  -  4(7^] 

(3^3  -  A*  +  4a^^„)* 


1600^2(0^  -  5„)tanh[Vi(Qm  -  0n)^]  -  4(7g] 
_ (3go^  -  A^„)* 

8<7o(v'2Am„tanh[\/2A 

mn  6]  -  4(7*1 
o\2  o  (®m  ^r»)  1 

*  ^  3-^  -  2A^„  -  3 - -j - } 

16  <7q  (7q 

The  three  equations  (2.8)-(2.9)  have  the  following  set  of  eight  solutions 
Pure  #1  ;  Am  -  1 1  A12  =0  a,  -  Irai 


Pure  4^2  :  /in  =  0  >ti2  =  ±1  <^2  —  “  (2.136) 

£ 

Mind  /111  r  .4i.>  -  r  aT  -  -  -(ai  *  0-2)  (2.13c) 

\  2  V  2  ‘  4 

Th<*  pure  mod*'.*^  r<)rresf>oiid  i<>  the  nuxlrs  of  the  .simple  eigenvalues  that  rualesce  to  form 
the  double  point.  They  share  the  same  natural  frequency  and  are  spatially  (horizontally) 
symmetric.  When  the  relevant  parameters  are  substituted  the  amplitude  correction  to  the 
natural  frequency  <7^  -  -  |oi  agrees  with  the  correction  found  for  the  simple  eigenvalues 
(equation  (2.5)  with  ^  =  1).  The  other  four  solutions  are  mixed  modes.  The  leading  terms  are 
proportional  to  the  sum  or  difference  of  the  two  pure  eigenfunctions.  The  amplitude  correction 
of  the  frequency  for  the  mixed  mode  solutions  <7™  differs  from  that  for  the  pure  modes  by 
an  amount  ^(oi  -  <12)  The  bifurcation  for  either  the  pure  or  mixed  modes  will  be  sub-  or 
super-critical  depending  on  the  value  of  S. 

Figures  la,b,c,d,  and  e  are  bifurcation  diagrams  for  the  natural  frequency  for  the  multiple 
eigenvalues  occuring  in  a  square  cross-section  with  mode  numbers  m  =  1  and  n  =  2.  The 
five  frames  show  the  effect  of  the  parameter  S  on  the  behavior  of  the  solutions.  Figure  la 
corresponds  to  infinite  depth  and  agrees  with  the  result  in  Bridges  (1986).  The  mixed  branch 
apparently  has  a  higher  natural  frequency  for  the  range  of  amplitudes  considered  and  for  all 
values  of  the  depth.  The  remainder  of  the  Figure  1  set  correspond  to  decreasing  values  of  6. 
As  S  decreases  the  branches  all  shift  to  the  right  and  eventually  (at  the  critical  depth  6*)  shift 
from  sub-  to  supercritical.  At  6  =  .115  there  is  the  interesting  property  that  the  pure  branch 
bifurcates  subcritically  and  the  mixed  branch  bifurcates  supercrilically. 

In  Bridges  (1986)  distributions  of  the  wave  height  for  the  pure  branches  and  the  mixed 
mode  branches  for  6  —  oc  are  shown  It  is  expected  that,  qualitatively,  the  distributions  for 
finite  S  (and  suitably  restricted  ()  will  be  similar. 

I  3.  Secondary  Bifurcation  when  t  <  0(6*) 

In  this  section  solutions  which  bifurcate  at  finite  amplitude  from  the  branches  emitted  by 
the  simple  eigenvalues  are  considered.  Therefore  a  perturbation  is  added  to  the  known  primary 
branch  solutions  ($,  T,a;)  found  bv  Verma  &:  Keller, 

I 

;  (3.1a) 


6 


Figure  I.  Ftfcel  of  6  on  the  hifurra- 
tioii  from  t)ic  double  eigenvalues  for 
(  =  I  and  (rn,»)  ranging  over  1,2. 
The  dashed  lines  correspond  to  the 
simple  modes  (left  branch)  and 
03,3  (right  branch)  and  the  solid  lines 
are  for  pure  (left  branch)  and  mixed 
(right  branch)  modes  emitted  by  the 
double  point  0|  3,  03,1. 
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These  expressions  are  substituted  into  the  governing  equations  and  boundary  conditions  which 
are  ihen  lineanxcd  about  (he  known  primary  branch  solutions.  The  linear  jirobleni  to  be  solved 
for  the  points,  if  thex  exist,  of  secondary  bifurcation  is 


dx^  ^  dy^  ^ 


and  on  y  = 


a<  di  dx  ^  di  di  dz  dx  ^  dz  dz  i  dy 


d0'  ,rd4>d<A'  d«‘d(^'  d4>d(/)'i  ,  „ 

and  in  addition  the  normal  derivative  is  required  to  vanish  on  the  solid  boundary.  This  is 
a  linear  differential  eigenvalue  problem  with  known  nonconstant  coefficients.  However,  the 
specific  value  of  f'  where  the  secondary  bifurcation  takes  place  is  sought.  Therefore  e'  is  the 
eigenvalue.  Since  c'  appears  nonlinearly  it  is  a  nonlinear  in  the  parameter  eigenvalue  problem 
and  there  is  the  further  complication  that  e'  is  responsible  for  the  size  of  the  domain.  The 
qualitative  shape  of  the  domain  is  known  since  T(i,z,<)  is  a  known  function,  but  the  precise 
amount  of  T(i.z,f)  is  the  unknown  eigenvalue.  This  is  to  be  constrasted  with  the  original 
eigenvalue  problem  in  Section  2  where  w  was  the  eigenvalue,  e*  was  a  variable  parameter,  and 
the  shape  of  the  free  surface  (and  hence  the  domain)  was  an  unknown  function. 

To  solve  this  eigenvalue  problem,  the  conjecture  of  Bauer,  Keller,  Sc  Reiss  (1975),  that  a 
secondary  bifurcation  may  occur  in  the  neighborhood  of  a  multiple  eigenvalue,  is  used.  For 
brevity,  the  analysis  is  undertaken  in  the  neighborhood  of  ^  =  1.  It  is  expected  that  a  similar 
analysis  will  hold  in  the  neighborhood  of  other  values  of  ^  at  which  double  eigenvalues  occur. 

It  was  shown  in  Section  2  that  for  ^  =  1  there  is  a  double  eigenvalue  for  every  pair  (m,n) 
such  that  m  3-  n.  At  the  double  eigenvalue  the  bifurcation  points  arc  given  by  (2.1)  with  (2.6). 
In  the  neighborhood  of  ^  -  i  this  double  eigenvalue  splits  into  two  primary  branches  emitted 
by  the  bifurcation  points 

n  =  f^n^  lanh  (3-5) 


(^n,m  ~  Y  Jr\/n*  4-  tanh[w6  4-  (3-6) 

A  measure  of  the  neighborhood  of  ^  -  1  is  given  b\  the  small  parameter  //  defined  by 

i  ■  1  -+  where  T  iitgn(^  -  1)  -  ::  1  (3-7) 

Following  the  conjecture  of  Bauer,  Keller,  ii  Reiss  that  the  secondary  bifurcation  disapears  at 
the  double  eigenvalue,  the  point  t'  on  the  primary  branches  where  the  secondary  bifurcation 
will  take  place  is  expressed  as 

('{n)  =  bon  ■*-  b\ti^  +  •■•  (3.8) 

The  solutions  previously  obtained  by  Verma  ie  Keller  for  the  primary  branches  emitted  by 
simple  eigenvalues  are  recast  using  (3.7)-{3.9)  as  expansions  in  the  small  parameter  fi, 


^  =  tpQ  +  + 


T  =  To  +  mT,  + 


(jJ  —  u^o 


(3  9a) 


{3  9b) 
(3.10) 


A  separate  analysis  is  undertaken  for  each  of  the  primary  branches  and  The  necessary 
details  for  the  analysis  along  the  branch  will  be  given  and  the  result  only  will  be  stated 
for  the  <r„„  branch. 

Substituting  the  expressions  (3.7)-(3.10)  into  the  linear  eigenvalue  problem  and  postulating 


<t>'  -  n<t>\  +  + 


r}'  -  +  ••• 


(3.11a) 

(3.116) 


Expanding  the  free  surface  boundary  conditions  in  a  Taylor  series,  and  equating  terms  propor¬ 
tional  to  like  powers  of  /i  to  zero  results  in  a  sequence  of  boundary  value  problems.  The  fact 
that  Wo  is  a  double  eigenvalue  results  in  the  leading  term  in  the  set 

w,,  ::  Afi  tanh  Ao6  (3  12) 


rii{Xy  2,t)  ~  All  COSamXCOS&nZ  4^  A\2C0S  l3f^XCOSan^Z\  sin  / 


(3.14) 


,  1  cosh  Am„(y (5) 

(i, y, 2, t)  = - — — — - - —  >lii  rosQ,„aToso„2  ~  /4)2Cos;?„icosam2  cos<  (3.15) 

u/’ti  cosh  A-mn^ 


and  the  normalization  is  taken  to  h 


■4tl  ^  Al,=:\ 


(3.16) 


The  problem  is  carried  in  the  usual  way  to  higher  order.  At  third  order  application  of  the 
double  solvability  condition  results  in  the  equations 


blAu  =  0 


(3.17) 

(3.18) 


which  with  (3.16)  form  a  set  of  three  equations  for  the  three  unknowns  bo.  An,  and  Ajj.  The 
term  03  is  given  by 

os(<)  -  J55(2?A,+  -  24^  -  9^)  +  8„J 

^  (3u>«  +  a;  -  4<,lf  + 

32w^(aTO  tanh  2a,„6  -  2u)l)  ^  32u;^(/3n  tanh  2/3„6  -  2w^) 

_ Ho  -  (ctm  -  ;3n)^l|4a;g  -  Ag  +  8Q^/?n] 

16w^{-4w^  +  Vi{0‘m  -  /5„)tanhl>/2(Q„i  -  /S„)6]} 

_ [i^O  -  (Qm  +  0„Y][4o}^  -  Ag  -  Botm0n\ 

16w*{-4u;^  +  \/2(a„,  +  /3„)  tanh[\/2(a„  +  /Sr.)^!} 

_ 1^0  ~  ~  -^O] _  jg. 

8wg{  — 4ooq  4-  v^Aq tanh(v/2Ao^]} 

For  sufficiently  large  b  it  htis  been  shown  by  numerical  evaluation  for  m,  n  ranging  over  1  to  10 
that  the  expression  for  03  is  positive  definite.  In  the  limit  as  6  — *  0  however  03  — *  -  jfg  Ag^ 
Therefore  as  6  —*  0  03  will  eventually  become  negative.  For  example  when  (m,n)  =  (1,2)  03 
changes  sign  when  6  ~  .075.  When  changes  sign  this  means  that  the  secondary  bifurcation 
will  “jump”  from  one  branch  to  another  which  would  be  a  discontinuous  phenomena.  It  is 
more  likely  that  this  behavior  is  a  ramification  of  the  non-uniformity  in  b  of  the  solution. 

The  solutions  to  (3.16)-(3.18)  are 


C ci8€  I  ;  6g  -  0,  .4 1 1  -  1,  .4)2  0 


Case  II  :  bn 


^ 1  .  T  .4 , .  =  n  A., 


(3,20a) 


1  onfct 


The  solution  in  (3.20a)  shows  that  the  basic  solution  bifurcates  from  the  primary  branch.  In 
(3.20b)  the  solution  for  the  secondary  bifurcation  on  branch  is  given.  The  ±  sign  shows 
that  the  bifurcation  takes  place  in  both  the  upper  and  lower  c  half-planes.  The  jump  to  the 
.4i;>  ^  0  solution  is  often  referred  to  as  mode  jumping  because  the  solution  acquired  on  the 
secondary  branch  is  qualitatively  different  from  that  on  the  primary  branch.  The  radical  in 
(3.20b),  when  as  >  0,  requires  that  (q„  -  ^„)r  <  0  for  secondary  bifurcation  to  occur  on 
branch  ffmn- 


A  similar  analysis  for  the  branch  results  in  the  bifurcation  equations 

--^(•5.  -  +  (AS  -  =  0 

Uln  ZAn 


bnAli  —  0 


(3.21a) 

(3.216) 


which  gives  the  points  of  secondary  bifurcation  on  that  branch.  For  convenience  define  Cm.n  to 
be  the  point  of  secondary  bifurcation  on  the  am„  branch  and  Cn.m  to  be  the  that  on  the  Onm 
branch,  then 

«m,n(M)  =  6m.nM  +  0(M*)  (3.22a) 

«n,m(M)  =  6n.mM  +  0(M*)  (3  226) 

Retaining  the  positive  branch  only  for  brevity  the  bifurcation  equations  on  each  branch  show 
that  _ 


6m, n  — 


6m, n  — 


2A*a3(6) 


)[l^(Ag  u;<) 

(3.23a) 

Wo 

|l^(Ag-^)^] 

(3.236) 

~  V  2A?as(6) 

A  secondary  bifurcation  takes  place  on  one,  and  only  one  at  a  time,  branch.  Noting  that 
T  =  signf^  -  1)  the  branch  on  which  the  secondary  bifurcation  takes  place  is 


{  -  1  (om  -  Branch 


(n,m) 


(m,n) 


(m,n) 


(fi,m) 


In  summary,  as  ^  departs  from  ^  =  1,  the  split  primary  bifurcation  points  given  by 
(3.5),  (3.6)  move  away  from  the  double  point.  When  $  >  1  they  both  move  to  the  right 
and  when  f  <  1  they  boili  move  to  the  left.  However  in  all  four  cases  given  in  the  table  the 
secondary  bifurcation  takes  place  on  the  branch  which  is  emitted,  after  splitting,  by  the  largest, 
in  magnitude,  of  the  two  bifurcation  points,  regardless  of  the  sign  of  r. 

It  has  been  shown  that  in  the  neighborhood  of  a  square  cross  section  ($  =  1)  a  secondary 
bifurcation  will  occur  on  one  of  the  two  branches  emitted  from  the  simple  eigenvalues  which 
result  from  the  splitting  of  the  double  eigenvalues.  By  expanding  in  the  neighborhood  of  this 
point  an  asymptotic  representation  of  the  solution  along  the  secondary  branch  may  be  found. 
A  small  parameter  i/  is  defined  as  a  measure  of  the  distance  from  the  point  of  secondary 
bifurcation.  It  is  assumed  that  the  parameters  are  such  that  the  secondary  bifurcation  point 
occurs  on  the  Cmn  branch.  A  similar  analysis  may  be  performed  for  a  bifurcation  from  the 
f7„„  branch. 

A  perturbation  is  added  to  the  known  primary  branch  solution 

4-  (3.24fl) 

ri  =  (r  +  Tj'  (3.246) 

w  —  ^  4“  Q  (3,24c) 

and  the  unknown  solutions  on  the  secondary  branch  are  expressed  as  a  regular  perturbation 
series  in  i/  and  p 


<i)'  ~  t'(p0ii  4-  p*<6i2  4-  •  •  •)  4-  t'*(p^2i  +  4-  •••)  +  "  ■  (3.25) 

T]'  -  I^(prj„  +  +  ■■■)  +  4-  •••)  +  ••  •  (3-26) 

n  —  t'(nio  4-  pHii  4-  4-  •  •  •)  4-  4"  4"  ^*n22  4"  •  •  •)  4"  •  ■  •  (3.27) 

The  substitution  of  expressions  (3.24)-(3.27)  into  the  governing  equations  and  boundary  con¬ 
ditions  results  in  a  set  of  boundary  value  problems  for  the  unknowns  <6,j  and  r],j.  The  analysis 
although  straightforward  is  lengthy  and  the  details  will  be  omitted. 


The  first  order  in  p  problem  results  in  =  0  for  all  j  and 

1  cosh  Ao(t/  4-6)  _ 

<f>u  ~ - ~rT~l — cos cos Q„r cost 

Wo  cosh  A(,6 


(3.28a) 
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»Jll  =  COS^nZCOSQm^sint  (3.286) 

and  the  higher  order  terms  (in  u)  are  omitted  for  brevity.  The  problem  of  order  i/*  results  in 
fl'id  -  ll-ji  —  0,  021  —  't-i\  '  and 


ll-jl  —  0,  021  -  't-i\ 

n-1'2  1 

-^0 


32 


IK 


'3w,i 


16aij^[4aiQ  -  \/2Antanh(v^2Afi6)' 


+ 


[3w3  -  A^  +  4am3ni  '‘ 


32w2l4w2  -  >/2(a„  -  ^„)tanh(v/2(Q„  -  3n)S)\ 


[3u;^  -  A^  -  Aoim3n 


32w^[4a;^  -  >/2(a„  +  ^„)tanh(%/2(a„,  +  /S„)6)j  ^  ^ 

and  the  other  higher  order  terms  are  omitted.  The  result  (3.29)  provides  an  expression  for  the 

frequency  along  the  secondary  branches.  The  complete  expression  for  the  natural  frequency  in 
the  neighborhood  of  the  double  eigenvalue  is 

<jJ  —  <jJq  “t"  4“  "i”  ^(/^^»^^)  (3.30) 


Since  H]]  is  proportional  to  quadratic  terms  the  sign  of  022  determines  whether  the  bifurcation 
is  sub-  or  supercritical.  Although  no  proof  has  been  undertaken  the  following  points  regarding 
the  stpn|n23j  are  made  based  on  numerical  evaluation  of  (3.29).  As  6  oo  the  st9n|n22j  <  0 
for  all  mode  numbers.  As  6  is  decreased  a  critical  value  of  6  is  reached  where  022  changes  sign 
and  this  critical  value  differs  for  different  mode  numbers.  For  example  when  (m,n)  =  (1,2) 
O22  changes  sign  from  -r  to  -  when  6  0.1155  resulting  in  a  shift  of  the  secondary  bifurcation 

from  sub-  to  supercritical.  Since  the  critical  value  of  6  on  the  primary  branch  is  slightly  different 
from  the  critical  value  for  the  secondary  branch  there  is  a  small  range  of  6  where  the  primary 
branch  is  supercritical  and  the  secondary  branch  is  subcritical. 


Figures  2a,b,c,d,  and  e  give  an  illustration  of  the  secondary  bifurcation  phenomena  for 
various  6  when  (m,n)  =  (1,2).  Figure  la  for  6  =  .20  is  similar  to  the  infinite  depth  result 
obtained  in  Bridges  (1986).  The  remainder  of  the  sequence  in  Figure  2  shows  the  shifting 
of  the  branches  to  the  right  as  6  is  decreased.  In  Figure  2d  the  secondary  bifurcation  is 
almost  vertical  as  022  0  here,  and  in  Figure  2e  the  secondary  bifurcation  has  shifted  to 
supercritical.  An  example  of  the  distribution  of  the  wave  height  (for  6  — *  00)  as  the  solution 
shifts  from  the  primary  to  the  secondary  pranch  is  shown  in  Bridges  (1986).  The  wave  held 
becomes  more  complex  as  the  solution  on  the  secondary  branch  is  acquiud.  The  addition  of 
the  hnite  depth  is  not  expected  to  signihcantiy  alter  qualitatively  this  distribution  for  suitably 
restricted  amplitude. 


■ 


m 
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Figure  2.  Effect  of  b  on  the  secondary 
bifurcation  when  (m,n)  =  (2, 1)  and  (  = 
|.  After  the  splitting  the  secondary  bi¬ 
furcation  occurs  on  the  oj.i  branch.  As 
b  is  decreased  the  primary  and  then  the 
secondary  bifurcation  shifts  from  sub-  to 
supercritical. 


2.7 


'VVV 


4.  Change  of  Type  When  <  =  0(6’) 

The  solutions  obtained  in  Sections  2  At  3  are  not  uniformily  valid  in  the  t  -  6  plane.  The 
higlter  order  terms  are  no  longer  of  higher  order  when  (  ~  'J’lierc  fore  in  this  section  a 

se  parate  analysis  is  performed  for  that  region  of  the  <  -  6  plane  where  f  -  0(6’).  This  is  done 
by  taking  6  =  <^/Tl  where  r  =  0(1)  and  carries  the  sign  of  t  (and  is  different  from  the  r  used 
in  Section  3).  The  governing  equations  are  (1.1)*(1.4)  but  with  the  scaling  modified  so  that 
6  appears  explicitly  in  the  equation  (y  is  scaled  with  h  instead  of  2a).  With  this  scaling  the 
linear  natural  frequency  has  a  finite  non-zero  limit  as  6  — >  0. 

Before  proceeding  to  the  fully  three  dimensional  problem  it  is  useful  to  analyze  weakly 
three  dimensional  waves  by  (following  Ablowitz  At  Segur  (1979))  considering  the  region  of 
parameter  space  where  {*  =  0(«),  or 

^  and  6  =  \/rt  (4.1) 

where  7  =  0(1)  and  carries  the  sign  of  e.  When  the  relations  (4.1)  are  substituted  into  the 
governing  equations  and  boundary  conditions  and  a  regular  expansion  in  t  is  sought  the  leading 
order  problem  is  a  wave  equation 

=  0  (4.2) 


where  U  is  the  D’Alembertian 


U  =  w2 


a’  a’ 


~  °at’  ax’  ' 

With  the  additional  requirement  that  vanish  on  the  vertical  boundaries  (4.2)  has  the 
general  solution 

Wo  =  am  (4.4) 

V'o(*,«,0  = /(f.*)  + /(X.^)  (4.5) 

where  Om  =  rnir,  f  t  X  =  ^  ~  am**  »nd  x  =  x  -b  The  leading  order  wave  height  is 

d  d 

7o(x,  z,  t)  =  -wo[^/(f,  z)  -b  ^/(x.  *)]  (4.6) 

the  unknown  function  /  is  found  through  application  of  the  solvability  condition  at  the  next 
order.  The  first  order  problem  is 

UV>,  =  F,(x,z,t)  (4.7) 


V  '#• 


r'vV 

,m  •  •  •  ■ 
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where  Fi  is  a  functional  of  zeroth  order  terms.  Solvability  of  (4.7)  requires  that 

r*2ck„. 


/ 


F, 


L:J(  ,  f  X 

2q,„  2 


rff  -  0 


(4.8) 


This  condition  is  derived  is  Bridges  (1985)  .A[)plicalion  of  (4.8)  results  in  a  (>artial  differential 
equation  for  /, 

u;o  a>:  ^  ^ 

When  -j  =  0  the  equation  for  /  can  be  integrated  to  yield 


fix)  ^  ^  +  Bcn^{x;K) 


(4.10) 


(where  A  and  B  are  constants)  the  usual  cnoidal  wave  and  subsequent  substitution  into  (4.5) 
results  in  a  standing  cnoidal  wave  (Bridges  (1985)).  With  the  retention  of  the  term 
the  equation  (4.10)  is  a  form  of  the  K-P  equation.  The  K-P  equation  was  first  derived  by 
Kadomtsev  id  Petviashvili  (1970)  in  their  study  of  the  stability  of  solitary  waves  to  transverse 
perturbations.  A  discussion  and  analysis  of  this  equation  can  be  found  in  Ablowitz  k  Segur 
(1979).  Dubrovin  (1981)  and  Segur  k  Finkel  (1984)  have  shown  that  this  equation  is  rich  in 
the  number  of  qualitatively  different  types  of  solutions  which  may  be  produced.  Here  the  right 
and  left  running  solutions  of  the  K-P  equation  would  be  combined  to  form  a  weakly  three 
dimensional  standing  wave.  One  example  of  a  solution  to  (4.9)  is  an  oblique  travelling  wave 
formed  from  f{x,  2)  =  /(x  “  0nz)- 

Instead  of  analyzing  this  equation  and  its  possibilities  further  an  analysis  with  (  unre¬ 
stricted  will  be  undertaken.  With  ^  unrestricted  and  6  =  y/rt  a  regular  perturbation  expansion 
in  (  is  assumed.  Substitution  into  the  governing  equations  and  boundary  conditions  results  in 
a  wave  equation  in  two  space  dimensions  at  leading  order. 

Wo'^-AV'o  =  0  (4.11) 


where 


dx^  ^  dz^ 


(4.12) 


and  it  is  required  that  the  normal  derivative  of  0o  vanish  at  the  vertical  boundaries.  The 
leading  order  wave  height  is  given  by 


r)o(x,z,t)  =  -wo^^V'o(*,2,0 


(4.13) 
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The  first  order  problem  results  in 


where 


u;'  --  )  --  V  ;/'i  ^  2,0 

ai-  o  ()t 


(4.14) 


buJi  dv’o 


For  solvability  it  is  required  that 


f^”  f  i  f  ‘i  dv/n 

III  ^~Gi(x,z,t)  dx  dz  dt  =  0 

Jo  J-iJ-i  ot 


(4.15) 


(4.16) 


The  function  ii)Q{x,  z,  t)  which  satisfies  this  functional  differential  equation  is  the  leading  order 
term  for  three  dimensional  standing  waves  in  a  rectangular  basin  when  e  =  O(S^).  A  general 
solution  to  (4.16)  has  not  been  found.  However  it  may  be  shown  by  substitution  that  with 


w.  =  -2ru;»;?^-  ^(2  -  ic»)i 


(4.17) 


one  poesible  expression  for  the  leading  order  wave  height  which  satisfies  (4.16)  is 


fJo(3t.2,0  =  h(t  +  Q„Z^d„z)^h(t  a„I-  -  Omi+/3„r)-f  fi(t  +  Q„Z-)?„2)  (4.18) 


where 


h(p)  -  A  -*■  Bcn*(p\K) 


(4.19) 


4  2,2  ,  E(jr;K) 

B  = 


(4.20) 

(4.21) 


where  E(x,k),  and  cn(/>;«)  are  Jacobian  elliptic  functions  (  Byrd  id  Friedman  (1971)).  Peri¬ 
odicity  in  time  and  the  finite  domain  require  that  k*,  the  modulus  of  the  elliptic  functions, 
satisfy  the  equation 

’-P  =  ° 

Jo  vl  -  K^sin  c 


Numerical  evaluation  results  in  '>■  0.9691. 


Therefore  a  solution  of  (4.16)  is  a  set  of  four  oblique  travelling  non-interacting  (to  leading 
order)  cnoidal  waves  which  when  combined  result  in  a  nonlinear  three  dimensional  standing 
cnoidal  wave.  It  is  also  illuminating  to  note  that  ih^  wave  height  (4.18)  may  be  expressed  as 


the  infinite  sum. 


r)o{x,z,t)  -  fco  cos cos 2  cos pt 


where  ko  is  a  constant  and 


.  _  1  [E(7r;Ac)  2  ,■] 

ao  —  — r - K  —  1 

/C*  L  JT  J 

2  pflP 

“p  ^  ^  t p  >  ° 
tc-*  1  - 


(4.23) 


(4.24a) 

(4.246) 


and  q  =  exp[-  K(  f ;  \/l  -  <c*)] . 


Examples  of  the  three  dimensional  cnoidal  standing  waves  are  given  in  Figures  3,4,  and  5. 
Figure  3  is  a  (m,  n)  =  (1,1)  mode  at  <  =  0,  Figure  4  is  a  (m,n)  =  (1,2)  mode  at  t  =  0,  and 
Figure  5  is  a  (m,n)  =  (2,2)  mode  at  t  =  0.  These  figures  give  a  prelude  to  the  richness  that  is 
possible  when  a  complete  solution  of  (4.16)  is  found. 

The  solution  obtained  in  (4.19)  is  only  one  possible  solution.  The  analysis  performed 
earlier  for  $  =  0(f)  which  resulted  in  a  K-P  type  equation  shows  that  at  least  in  the  region 
^  —  0(e)  there  are  solutions  other  than  the  four  travelling  oblique  cnoidal  waves  and  it  is 
expected  that  other  regions  of  f  will  also  have  additional  solutions. 
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additional  expansion  in  the  region  of  non-uniformity  shows  that  the  wave  field 
changes  type*  One  possibility  in  this  region  is  a  field  of  three-dimensional 
cnoidal  standing  waves. 


